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Abstract. The purpose of this brief paper is to prove De Giorgi type results 
for stable solutions of the following nonlocal system of integral equations in 
two dimensions 

L{ui) = Hi(u) in R 2 , 

where u = ( Ui )£h 1 for Ui : R n —»• R, H = (Hi)™ :1 is a general nonlinearity. 

The operator L is given by 

L(ui(x)) := / [ui(x) — Ui(z)]K(z — x)dz, 

Jr 2 

for some kernel K. The idea is to apply a linear Liouville theorem for the 
quotient of partial derivatives, just like in the proof of the classical De Giorgi’s 
conjecture in lower dimensions. Since there is no Caffarelli-Silvestre local ex¬ 
tension problem associated to the above operator, we deal with this problem 
directly via certain integral estimates. 
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1. Introduction 

We study the following semilinear elliptic system of nonlocal equations 

(1) L{ui) = Hi{u) in R n 

where u = (ui) r fE 1 for m : R n -A R and H = (Hi) 7 fL 1 is a sequence of locally 
Lipschitz functions. The operator L is an integral operator of convolution type 

(2) L(ui(x)) = / [ui(x) - Ui(z)]K(z - x)dz , 

J R n 

where the kernel K satisfies a certain number of conditions. We are interested in 
De Giorgi type results for this system in two dimensions with a general nonlinearity 
H. The case of scalar equations, that is when m = 1, is considered very recently 
by Ros-Oton and Sire in 24] and by Hamel and Valdinoci in [22j. 

Regarding semilinear equations, Ennio De Giorgi (1978) in [14j proposed the 
following conjecture that has been a great inspiration for many authors in the field. 

Conjecture 1. Suppose that u is an entire solution of the Allen-Cahn equation 
that is 


( 3 ) 


A u + u — u 3 = 0 


on 
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satisfying |w(x)| < 1, -§fh{x) > 0 for x = (x',x n ) G R”. Then, at least in dimen¬ 
sions n < 8 the level sets of u must be hyperplanes, i. e. there exists g G C 2 (R) such 
that u{x) = g{ax' — x n ), for some fixed a G R" _1 . 

Ghoussoub and Gui in [20] provided the first comprehensive affirmative answer 
on the De Giorgi conjecture in two dimensions. In fact their proof if valid for 
a general Lipschitz nonlinearity and not necessarily double-well potentials. Their 
proof uses the following, by now standard, linear Liouville type theorem for elliptic 
equations in the divergence form, which (only) holds in dimensions one and two, 
see ism If (f> > 0, then any solution a of 

(4) div(0 2 Vu) = 0, 


such that cjxj is bounded, is necessarily constant. This Liouville theorem is then 
applied to the ratio a := -Sr-/4^ to finish the proof in two dimensions. Ambrosio 
and Cabre in [J], and later with Alberti in (2| for a general Lipschitz nonlinearity, 
provided a proof for the conjecture in three dimensions by noting that for the linear 
Liouville theorem to hold, it suffices that 


(5) [ </>V < CR 2 , 

J B2r\Br 

for any R > 1. Then by showing that any solution u such that d Xn u > 0 satisfies 
the energy estimate 


( 6 ) 


f-a 2 < CR n 


Ib r 


they finished the proof. Even though the original conjecture remains open in di¬ 
mensions 4 < n < 8, there are various partial, yet groundbreaking, results for 
dimensions 4 < n < 8. In this regard, Ghoussoub and Gui showed in 21] that 
the conjecture holds in four and five dimensions for solutions that satisfy certain 
antisymmetry conditions, and Savin in [25| established its validity for 4 < n < 8 
under the following additional natural hypothesis on the solution, 


(7) lim u(x', x n ) —» ±1. 

tc n —>-±oo 

In [26| , Wang provided a new proof for the proof of Savin that involves more 
variational methods in the spirit. Unlike the above proofs in dimensions n < 5, 
the proof of Savin is non-variational and does not use a Liouville type theorem. 
We refer interested readers to sn for more information regrading this Liouville 
theorem for local semilinear equations. For the case of n > 9, del Pino, Kowalczyk 
and Wei in [l5[ provided an example that shows the eight dimensions in the critical 
dimension for the conjecture. In addition, for the case of system of semilinear 
equations, Ghoussoub and the author in [l8| established a counterpart of the above 
linear Liouville theorem for H -monotone solutions. See [16] for a slightly improved 
version of this theorem as well. Note that since u = {ui) 1 fL l is a multi-component 
solution, the concept of monotonicity needs to be adjusted accordingly. We borrow 
the following definition from [l8[ . 


Definition 1. A solution u = (uk)™=i °f I HP is said to be H-monotone if the 
following hold, 

(i) For every 1 < i < m, each Ui is strictly monotone in the x n -variable (i.e., 
d Xn Ui ± Oj. 
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(ii) For all i < j, we have 

( 8 ) djH i (u)d n Ui(x)d n Uj(x) > 0 for all x £ R“. 


Note also that (ii) implies a combinatorial assumption on the sign of partial 
derivatives of Hi and we call any system that admits such an assumption as ori- 
entable system. In other words, for orientable systems there exists a sequence of 
sign functions 6 = where 0* £ {—1,1} such that djHi{u)6i0j > 0. Coun¬ 

terparts of the above Liouville theorem for nonlocal equations with the fractional 
Laplacian operator are provided by Cabre and Sola-Morales in (Hi and Cabre and 


10, 111. In addition, for the case of system of equations Sire and the author 
in |19| and the author in [H] provided such a Liouville theorem. We refer interested 
reader to {§] by Cabre and Cinti for the fractional Laplacian case as well. 

This paper, in the light of the above results, deals with one-dimensional symme¬ 
try results for bounded stable solutions of m via proving a linear Liouville theorem 
for the quotient of partial derivatives. Note that for the case of fractional Laplacian, 
see [9h 12L |l9( and references therein, the extension problem that is a local problem, 
derived by Caffarelli and Silvestre in [ll|, plays a key role to prove one-dimensional 
symmetry results. But, such an extension problem is not known for a general op¬ 
erator L given in {2|). Therefore, we apply the integral definition of the operator L 
directly in our proofs. Most of the ideas and computational techniques applied in 
this work, in regards to the operator L, are strongly motivated by the ones applied 
in fH, [iij • Here is the notion of stability. 


Definition 2. A solution u = of {!]) is called stable when there exist a 

sequence of functions (j) = and a constant A > 0 such that each (f>i does not 

change sign. In addition, (j> satisfies the following 

m 

(9) L(cf>i) =’^2 / djHi(u)(j>j + \(fi in M™, 

j =i 

where djHi(u)<f)j(j)i > 0 for all i < j when 1 < i,j < m. 


It is straightforward to see that any U-monotone solution is a stable solution via 
differentiating H|) with respect to x n and defining = d Xn Ui. As it is shown in [H, 
in the absence of ff-monotonicity and stability there are in fact two dimensional 
solutions for Allen-Cahn systems in two dimensions. See also Q by Alikakos for 
a comprehensive note regarding semilinear local systems. Therefore, the concepts 
of LZ-monotonicity and stability seem to be crucial in the context. Let us fix the 
following notation. 


Notation 1. Throughout this paper we refer to sequences of functions if = (ipi) 7 fL 1 
and o = for ipi := V?q • r) where rj(x) = rj(x', 0) : R n_1 —»• R and <Tj := ^ 7 . 

Note that there is a geometric approach regrading the De Giorgi’s conjecture in 
two dimensions, but in this paper our technique is applying Liouville theorems. The 
next definition is the notion of symmetric systems, introduced in [H|- The concept 
of symmetric systems seems to be crucial for providing De Giorgi type results for 
system m with a general nonlinearity. 

Definition 3. We call system m symmetric if the matrix of partial derivatives of 
all components of H = ( y Hi) r fL l that is 

H := 
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is symmetric. 

Occasionally, we assume that the kernel K satisfies the following conditions as 
this is the case in j24j. 

(A) The kernel K(y ) is nonnegative, even and measurable function. In ad¬ 
dition, K has compact support in B\ satisfying K = 0 in R” \ B\ and 

f Bl \y\ 2K (y) d y < c - 

(B) Harnack inequality: Suppose that / does not change sign and is a solution 
of L(f) + g(x)f = 0 in R” for some g 6 L°°(R n ) then sup Bi ( Xo ) / < 
C inf Bl ( Xo ) / for any x 0 G R”. 

(C) Holder estimate: suppose that / is a bounded solution of L(f) = g in R” 
where g is bounded in B\. Then ||/||c«(B 1/2 ) <C(\\g\\ L ~ Bl +\\f\\L~mn), 
for some positive constants a, C. 

Note that the classical De Giorgi’s conjecture with the stability assumption, 
instead of monotonicity, is known as the stability conjecture and it is known to be 
true for n = 2 and not to be true in dimension n = 8. For two dimensions the 
same proof for monotone solutions work for stability conjecture and for the eight 
dimensions we refer to [23j by Pacard and Wei. This conjecture remains open for 
other dimensions. As a matter of fact our results in this work answer the stability 
conjecture in two dimensions for system of equations. 

2. Main results and proofs 

We start this section by deriving an integral system of equations for the sequence 
of functions cr = (af) J fL 1 where each eq is the quotient of partial derivatives, see 
Notation [l] 

Lemma 1. Suppose that u = ® s a stable solution of m■ Then for each 

i = 1, • • • , m and for any x G R” we have 

,, m 

(10) / [a t (x) - <Ji(z)\<t>i(z)K(x - z)dz = ^ d j H l (u)(j)j(x)[aj(x) - <Ji{x)\. 

Jr " 1=1 

Proof. Since u is a stable solution, there exist a sequence of functions </> = 
and a constant A > 0 such that 

m 

(n) £(<Mz)) = E djHi(u)(f>j + A cji- 

i=i 

Differentiating CD we get 

m 

(12) M^)) = E djHi(u)ipj + A ifi. 

l=i 

Since ifi = from m we have 

m 

(13) L(ai(x)((f>i{x)) ='^2d j H i (u)aj(x)(j)j(x). 

l=i 

Multiply HID with —<Ji and add with m to get 

m 

(14) L(ai(x)((j)i(x)) - a i L((j) i (x)) = ^ (x) - <Ji(x)). 

1=1 
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From the definition of the operator L we get 

L(<Ji(x)(4>i(x)) — tJiL(<f>i(x)) = [cri(x) — <Ti(z)]y>i(x) — <t>i(z)\K(a 

Jr 2 

= <f>i(x) / [<Ti(x) - cri(z)]K(x — z)dz 
J R» 

- / [<jj( x) - <Ji(z)][(t>i(x) - 4>i(z)\K(x - z)dz 

JR 2 

= / [<Ti(x) - <Ji(z)\(j)i(z)K(x - z)dz. 

Jr 2 


This and m gives the desired result. 


□ 


As an immediate consequence of Lemma [T] we provide the following integral 
estimate. 

Corollary 1. Suppose that assumptions of Lemma Q] hold. Assume also that rj = 
where r]i 6 C^(R n ) for each 1 < i <m. Then, 

m p p 

(15) ^2 / a i( x )[ a i( x ) - o-i(z)]<t>i(z)(t>i(x)K(x - z)t]f(x)dzdx 

i=1 JR n JR n 

m „ 

(16) = V' / djHi(u)<f>i(x)<f)j(x)[aj(x) - ai(x)\ai(x)r]?(x)dx. 

i,i= 

Proof. Multiply (fTUl) with <j>i(x)r]f (x) and integrate with respect to x. Then take 
sum on the index i. This finishes the proof. 

□ 


Here is our main result. 


Theorem 1. Suppose that u = (ui) r [L l is a bounded stable solution of symmetric 
system m in two dimensions where the kernel K satisfies conditions (A)-(C). Then, 
each Ui must be a one-dimensional function for i = 1, • • • , m. 


Proof. From Corollary [T| we have 

m ,, 

V / djHi{u)4>i{x)4>j(x)\(jj{x) - <Ji(x)\<Ji(x)rjf(x)dx 
ij=l jRn 

171 p p 

= V] / / cn(x)[<Ji(x) - cri(z)](f>i(z)(j>i(x)K(x - z)r]i(x)dzdx 

i-1 

171 p p 

= -^2 / ^i{z)[<Ti(x) - ai(z)]4>i(x)(l)i(z)K(x - z)r]i(z)dzdx. 

i=1 JR n J R n 


This implies that 


171 p p 

I = y2 / / \rft(x)cji(x)-rft(z)<Ji{z)][cri{x) - (Ji{z)\(t)i{x)(l)i(z)K(x - z)dzdx 

l=1 JR n J R n 

m ,, 

= 2 ^ / djHi(u)<j>i(x)(j>j(x)[<jj(x) - <Ji{x)\ai(x)rjl(x)dx. 

ij=1 jRn 


— z)dz 
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Set rji = (r where (r £ C“(K") is the standard test function that is £r = 1 in 
Br and Cr = 0 in \ B 2 R with ||VCr||l°“(b 2 b\-Br) < CR -1 . Note that for stable 
symmetric systems we have 


- Oi) = '^2,MjdjH i {u)cTi{(Tj - <Ji) + y (j)i4>jdjHi{u)(Ji{aj - cq) 
i,j i<j i>j 

= djHi(u)cri{Uj - <T») + y (j)j(j)jdjHi{u)aj(ai - aj) 


i<j 


i<j 


= ~y <j > i < t > jdjHi(u)(<Tj - CTi) 2 < 0. 

i<j 

Therefore, I < 0. Note that 

(!7) 2 [r)t(x)<Ji{x)-rjf(z)<Ji(z)\ = [<Ti(x) - cn(z)][r) 2 (x) + r] 2 (z)] 

(18) +[cri( x) + a, ( 2 )] [rj 2 (x) - rj 2 (z)}. 

From this and the fact that I < 0, we get 

J(R) [ai(x) - (?i(z)] 2 [(R{x) + C%(z)]<f>i(x)<j>i(z)K(x - z)dzdx 

&i(z) - <?i (*)][Cfl(®) - Cr(z)]M x )Mz)K{x - z)dzdx 

[<Ti(z) - (Ji(x)][(Ji(z) + (Ji{x)\ 

[C«(®) - (r{z)][(r{x) + C R(z)\<l>i(x)<l>i(z)K(x - z)dzdx. 

Note that J(R ) > 0. From Cauchy-Schwarz inequality, we have 



J\R) 



*) 




[<yi{x) - < 7 i(z)] 2 [(R{x) + CR{z)] 4 >i{x)(j)i{z)K(x - z)dzdx 

[cri(x) + ai(z)] 2 [^R(x) - C,R(z)] 2 (j)i(x)<j)i(z)K(x - z)dzdx \ . 

From the definition of (r it is straightforward to see that 

(C r( x ) - Cr(*0) 2 < CR~ 2 \x - z\ 2 Xb 3R (x) when \x - z\ < 1. 

Note also that from Condition (A) we have K = 0 on K 2 \ B\. Therefore, 

m r. r. 

( 19 ) y] / / [vi{x) + ai(z)] 2 [C,R{x)-C,R{z)) 2 (t>i{x)(t>i(z)K(x- z)dzdx 

, , JRri J-^n 


(20) < R 


1=1 

-2 


[<Ji{x) + (Ji{z)Y(t>i{x)(j)i{z)\x — z\ 2 K(x — z)dxdz. 


I —1 J B 2 R J B\ 


Applying the condition (C) to each equation in (JTJ we get |Vztj| £ L°°(R n ) for each 
i. From the definition of cq then we have |<7j| < ||V'Ui||£oo( R ,»)|(/> i | _1 pointwise in 


l . Therefore, for each * = !,••• , to, 


[<Ti(x) + di(z )] 2 < ||Vlti|| L oo( R n) (\(/>i(x)\ 2 + \(t>i{z)\ 2 ) . 
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This implies that the function in the integrand of (l20l) is bounded by 

MaO + ai(z)] 2 Mx)Mz) < IIVuillz^Rn) ’ 

where the latter in bounded due to the condition (C). This and CHI) implies that 

( 21 ) 

( 22 ) < 

(23) < 


in p p 

/ l 0 *^) + <Ti{z)\ 2 (j)i{x)(l)i(z)\x - z\ 2 K{x - z)dxdz 

l—l J B2R J Bi 

m p 

CR-zy' / / \x — z\ 2 K(x — z)dxdz 

i = l J B2R J Bi 


CR~ 2 \B 2R \, 


where we have used the condition (A) in the last inequality. Note that in two 
dimensions R^ 2 \B 2 r\ is bounded uniformly in R. This implies that the integral in 
(fl 9 l) is bounded uniformly in R, i.e., 



W i(x) + cr i (-2)] 2 [Cfl(a;) - (R(z)] 2 (/>i(x)(/>i(z)K (x - z)dzdx < C , 


where C is independent from R. This implies that J(R) < C for any R and therefore 

m p p 

V / / [<Ti(x) - ai(z)] 2 (/)i(x)(l)i(z)K(x - z)dzdx < C. 

i=1 JR n JR n 

The above upper integral bound on J shows that 

m p p 

^ - (Ji{z)] 2 (j)i{x)(j)i{z)K{x - z)dzdx = 0. 

;_i J R n JR 71 


i =1 


The fact that each </>, does not change sign guarantees that each eq must be constant. 
From the definition of cr; we get Vu; • rj = Ci<pi in R 2 . Therefore, each Ui is a one¬ 
dimensional function. 

□ 


For the rest of this section, we analyze the stability of solutions. For the semi- 
linear equation — Au = f(u) in R", it is known that the pointwise stability that 
refers to the linearization of the operator implies a stability inequality of the form 

(24) / f(u) C 2 < / |VC| 2 , 

Ji” 

for any test function £ £ C^R”). This inequality has been extensively used in the 
literature to study qualitative behaviour of solutionsjrf the above equation. As a 
matter of fact, it is shown by Ghoussoub and Gui in [20], see also [l] by Berestycki, 
Caffarelli and Nirenberg, that the stability inequality also implies the pointwise 
stability. This equivalence property is known for the case of system of semilinear 
equation by Ghoussoub and the author in [Til as Lemma 1. In addition, for the 
case of integral equations of the form ©, when m = 1, such equivalence is also 
known by Ros-Oton and Sire in [24}. In the next lemmata we provide a similar 
result for system CD- 

Here is the stability inequality that is a counterpart of (l24l) . 
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Lemma 2. Let u = denote a stable solution of symmetric system fTp. Then 

771 n -i Tin n rt 

(25) "52 / djHi(u)Ci(x)Cj(x)dx < - V] / / [Ci(z) - CiOO] 2 -^ - x)dzdx, 

; ;_i JR" ^ • , JR" JR" 


*|J = 1 


/or any £ = (COT where Q € C/°(K n ) for 1 < i < m. 

Proof. Suppose that u = (ui)^L 1 denote a stable solution of HJ. Then there exist a 
sequence of functions (j> = (</fc)fcLi and a constant A > 0 such that djHi(u)<f>j<f>i > 0 
and 


(26) 


u)(j)j + A cf>i in 
i=i 


c 2 

Multiply both sides with where / = (Ci)£Li is a sequence of test functions. 
Therefore, 

c 2 

V Si 


(27) 


E L (^)x = E Wu)xC?+ac? 

2,.7=1 


2 = 1 


Note that the left-hand side can be rewritten as 

^2 m z-2 n 


i,. 7=1 


i<j 


i>3 


E djHitoh j = x)w«)^x + E^(“)^^+Ew u K' 

2 \ rn 

L + E diHi ( U )ti 

i J 2=1 

> 2 £ i///,(u)QCi + E W«Ki 


m / d 2 C: 

= E ( ^ H ' ( M ) ^ X + ^ 


i<j 


i<j 


2=1 


= E 'V/'Oogcj- 

»d=i 

From this and (HZD we get 


771 r. 771 « 

E / djHi(u)CiCj < E / L (^) 

, =i Jr« i=1 7 r« 


(28) 

bi=i 

Note that for each i we have 

(29) f L((f) i )^- = ^f j [<j)i(x) -<fi{z)) 

J R n Vi z JR n ./R n 


c 2 

Si 


C 2 (^) C?(z) 


K(x — z)dxdz. 


/R" JR” 

The rest of the proof deals with simplifying the integrand of the right-hand side. 
To do so, we closely follow the computational techniques provided in 23,124], Note 
that 

’C 2 (*) C 2 W 


(30) [</>i(a;) - 0 i(z)] 


_4>i{x) (fi{z)_ 


= 2[4>i{x) - (j>i(z)\ [Ci(a;) - Ci(~)] 


C i(x) + (i(z) \ ( [4>i(x) + (j)i{z)]‘ 


4>i(x) + 4>i(z) J \ 4(f>i(x)(f>i(z ) 


~[Mx) - Mz)? 


Ci(x) + (i(z)\ f 2(i(x) + 2(f(z)\ f[(f>i(x) + M Z )\ 


<t>i{x) + <t>i{z) J V [Ci(*) + (i( z )] 2 J V 4 <j>i(x)<t>i(z) 
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The fact that for stable solutions, for each i function (pi does not change sign implies 
that 4>i(x)(f>i(z) = \4>i{x)(j)i{z)\. On the other hand, we have 


2 [M x ) ~ <&(*)] [Ci(*) - <i(z)] 


Ci(*) + C»fc) 

<f>i{x) + (j)i(z) 


< [C»(*) - + [<t>i(x) - Mz)]- 


<,h(x ) + 4>i{z) 


From this, (ED and ED we get 

d 2 1 m 
,c 


/ R' 


This and 


L{<t>i ) 

^ z i=i ■ 

finishes the proof. 


[Q(x) - C,i{z)\ 2 K(z - x)dzdx. 


□ 


We now provide the other side of the argument to conclude the equivalence of the 
linearized stability (OBI) and the stability inequality (l?5l) . Note that the follo wing 
proof is strongly motivated by the one given by Ghoussoub and the author in 18] 
for the case of semilinear local systems. 

Lemma 3. Suppose that u = (ui) r fL 1 is a solution of symmetric system W for 
which the inequality if 25\) holds. Assume that conditions (B) and (C) hold as well. 
Then u is a stable solution of W- 

Proof. For each R > 0 define 
(31) 

A i(R) := inf J T R ( 0 ; < = (QZi where & G H K (R n ), & = 0 in R" \ B R , V f g 

{ i= i Jb R 

where Hk{B r ) is the closure of C“(K") with the norm 

ll/lk*( R ") : = / / I f(x)~ f(z)\ 2 K(x-z)dxdz, 

JR n JR n 

and the operator Tr is given by 

m 

(32) T r ( C) := Y] / d j Hi(u)C i (x)C j (x)dx 

ij=i jRn 
I m r r 

(33) / [Ci(x)-Ci(z)] 2 K(z-x)dzdx. 

Z i=1 JR n J R” 

From (l25l) we have \\{R) > 0 and there exist eigenfunctions (7 such that for each 

i = 1, • • ■ , m, 


(34) 


L(cn = Z7=1 W«)Cf + MW, if 1*1 < R, 

cf = 0, if jxj = R. 


Since the system is orientable, there exists (0 i ) r fL 1 such that djHi(u)6idj > 0. We 
now apply the maximum principle and use the sign of each 0i to assign signs for 
the eigenfunctions (£ 7)k so that they satisfy 
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Here £/* can be replaced by sgn{6i) | if needed and we can normalize each Q so 
that 

m 

(36) ElC^(0)l = l- 

i=1 

Note that Xi(R) f Ai > 0 as R —> oo. Define xf '■= s 9 n (Cf')Cp and multiply system 
(l34ll with sgn(Cf) to get that 

I = d i Hi ( u )x? + \ d o H i{u)\xf + X\{R)xf > if M < R, 

\ X? = 0, if \x\ = R. 

Here we have used the fact that from orientability, sgn((f- £f) = sgn(djHi(u )) 
when i ^ j. Each x f is a nonnegative function and x — {Xi)iL i is a solution for 
(EH). From the Harnack’s inequality, given in the condition (C), for any compact 
subset H, maxf! \xf \ < C'(fi) minn |xf| for alii = 1, • • ■ , m with the latter constant 
being independent of xf'- From elliptic estimates, given as the condition (B), one 
can see that the family (xf )r have uniformly bounded derivatives on compact sets. 
So, for a subsequence ( Rk)k going to infinity, (xf fc )fe converges in Cf oc (lR") to some 
Xi G C 2 (IR") and that Xi > 0. From (1371) we see that xi satisfies 

(38) L(xi) = d i H i {u)xi + 'Y^\djH i (u)\xj + \iXi 

> (diHi(u) + Ai)xi 

Note that Xi ^ 0 an£ f diHi(u) is bounded. The strong maximum principle implies 
that either Xi = 0 or Xi > 0 in R". Suppose that Xi = then from (1551) we have 
l^i^*( M )IXi = 0 that is Xj = 0 if j ^ i. This is in contradiction with (1551) . 
Therefore, Xi > 0 for all * = 1, • • • , m. Finally, set 4>i := sgn(9i)Xi for i = 1, ■ ■ • , m 
and observe that (4>i)i satisfies (1551) and that djHi{u)(j)j(j)i > 0 for i < j, which 
means that u = (ni)^Ti is a stable solution. 

□ 
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